Spontaneous synchronization of metronomes on a freely rotating platform 
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Spontaneous synchronization of an ensemble of metronomes placed on a freely rotating platform is 
studied experimentally and by computer simulations. The experiments revealed that independently 
of the initial conditions and metronome numbers on the platform, only in-phase synchronization 
is observable. The conditions favoring the emergence of this collective behavior is investigated. It 
is found that there is a limiting natural frequency of the metronomes below which spontaneous 
synchronization is not observable. Increasing the number of metronomes on the platform leads to a 
decreasing trend in the synchronization level and a longer time to achieve this. 



I. INTRODUCTION 

Spontaneous synchronization of coupled non-identical 
oscillators is a well-known form of collective behavior [l|- 
0]. The problem has been intensively studied since Huy- 
gens, and if the legend is true he was the first one who 
noticed and reported the synchronized swinging of pen- 
dulum clocks suspended on the same wall. By simple 
experiments he found that the synchronized state is a 
stable limit cycle of the system, because even after per- 
turbing the system, the pendulums came back to this dy- 
namic state. Originally Huygens thought, that this odd 
kind of sympathy, how he named it, occurs due to shared 
air currents between the pendulums. He performed sev- 
eral tests to confirm this idea. His experimental setup 
was really simple, with two pendulum clocks hanged on 
a common suspension beam which was placed between 
two chairs 0] ■ After performing some additional tests 
Huygens observed a stable and reproducible anti-phase 
synchronization, and attributed this to the impercepti- 
ble vibrations in the suspension beam. He summarized 
his observations in a letter to the Royal Society of London 
, and launched the study of synchronization phenom- 
ena and coupled oscillators. 

Recent studies aimed to reconsider in various forms 
Huygens' two pendulum-clock experiment and realisti- 
cally modeled the system. Bennet and his group Q in- 
vestigated the same two pendulum-clock system as Huy- 
gens did and came to the conclusion that as a function 
of the system's parameters several types of collective dy- 
namics is observable. For strong coupling usually a "beat 
death" phenomena appears, where one pendulum oscil- 
lates and the other does not. For weak coupling, synchro- 
nization does not appear, and a quasi-periodic oscillation 
is observable. There is however an intermediate coupling 
strength interval where the anti-phase synchronization 
observed by Huygens appears. Huygens had thus some 
luck with his setup, because the coupling was just in the 
right interval: strong enough to cause synchronization, 
but also weak enough to avoid the "beat death" phe- 
nomenon. 

Dilao [13] came to the conclusion that the periods of 
two synchronized nonlinear oscillators (pendulum clocks) 



differ from the natural frequencies of the oscillators. Ku- 
mon and his group [llj studied a similar system consist- 
ing of two pendulums, one of them having an exciting 
mechanism, and the two pendulums were coupled by a 
weak spring. Fradkov and Andrievsky (l^ elaborated a 
model for such a system, and obtained from numerical 
solutions that both in- and anti-phase synchronizations 
are possible depending on the initial conditions. 

Kapitaniak and his group [l^ revisited Huygens' orig- 
inal experiment and found that the anti-phase synchro- 
nization usually emerges, although in rare cases in-phase 
synchronization is also possible. They also elaborated a 
more realistic model for this experiment [l^ . 

Pantaleone [l^ considered a similar system, but he 
used metronomes placed on a freely moving base instead 
of pendulum clocks. He sir nply modeled the metronomes 
as van der Pol oscillators [la] and came to the conclu- 
sion that anti-phase synchronization occurs in some rare 
cases only. He proposed this setup as an easy classroom 
demonstration for the Kuramoto model and extended the 
study for larger systems containing up to seven globally 
coupled metronomes. Recently Ulrichs and his group [13 
examined the case when the number of metronomes was 
even larger. 

The present work intends to continue this line of stud- 
ies considering an ensemble of metronomes arranged sym- 
metrically on the perimeter of a freely rotating disk, as 
illustrated in Figure [T] The free rotation of the disk 
acts as a coupling mechanism between the metronomes 
and for high enough ticking frequencies synchronization 
emerges. Our aim here is to investigate the conditions 
favoring such spontaneous synchronization in large en- 
sembles of metronomes. In order to achieve our task, we 
first study the dynamics of the system by well-controlled 
experiments. Contrary with earlier studies that investi- 
gated the final stable dynamic state of the system, here 
we also consider and describe the transient dynamics 
leading to synchronization. We construct then a realistic 
model for the system and its modeling power is proved 
by comparing its results with the experimental ones. Fi- 
nally, the model is used to investigate synchronization 
level and emergence of synchronization in large ensem- 
bles composed of identical metronomes. 



2 



II. EXPERIMENTAL SETUP 

The experimental setup is sketched in Figure [T] The 
main units are the metronomes (Figure [T] and [Hh-), which 
are devices that produce regular, metrical beats. They 
were patented by Johann Maclzel in 1815 as a timekeep- 
ing tool for musicians ((3). The oscillating element of 
the metronome is a physical pendulum, which consists of 
a rod with two weights on it (Figure : a fixed one at 
the lower end of the rod denoted by Wi and a movable 
one W2, attached to the upper part of the rod. In general 
Wi > W2, and the rod is suspended on a horizontal axes 
between the two weights in a stable manner, so that the 
center of mass lies below the suspension axes. 




FIG. 1: Schematic view of the experimental setup: 
metronomes are placed on the perimeter of a disc that 
can rotate around a vertical axes. 



By sliding the W2 weight on the rod, the oscillation 
frequency can be tuned. There are several marked places 
on the rod where the W2 weight has a stable position, 
yielding standard ticking frequencies for the metronome. 
These frequencies are marked on the metronome in the 
units: Beats Per Minute (BPM). 

Another key part of the metronomes is the exciting 
mechanism, which has the aim to compensate the en- 
ergy lost by frictions. This mechanism yields additional 
momentum to the physical pendulum in form of pulses 
delivered at a given phase of the oscillation period. For 
a more detailed analysis of this exciting mechanism we 
recommend the work of Kapitanik ct. al [l^ . 

For the experiments we used the commercially avail- 
able Thomann 330 metronomes (Figure |3K) . From the 
10 metronomes we had bought, 7 ones were carefully se- 
lected, so that their standard ticking frequencies be as 
similar as possible. Naturally, since there are no two 
identical units, in experiments we have to deal with a 
non-zero standard deviation of their natural frequencies. 

In order to globally couple the metronomes, we placed 




FIG. 2: Schematic view of the metronomes' bob. The 
dotted line denotes the horizontal suspension axes, and 
the white dot illustrates the center of mass. 

them on a disk shaped platform which can rotate with a 
very small friction around a vertical axes, as it is sketched 
on Figure [1] and illustrated on the photo in Figure [3^. 

In order to monitor separately the dynamics of 
all metronomes, photo-cell detectors (Figure [SId) were 
mounted on them. These detectors are commercial ones 
(Kingbright KTIR 0611 S), and contain a Light Emit- 
ting Diode and a photo-transistor. They were mounted 
on the bottom of the metronomes. 

The wires starting from each metronome (seen in fig- 
ure [3l) connect the photo-cells with a circuit board, al- 
lowing data collection through the USB port of a com- 
puter. Data collection was done using a free, open-source 
program, called MINICOM. {^). The data is saved 
in log files, and can be processed in real-time. It is 
possible to follow simultaneously the states of up to 8 
metronomes. The circuit board sends data just in those 
moments when there is a change in the signal from the 
photo-cell system (i.e. a metronome's bob passes the 
light-gate). In such moments it records a string like: 
0-0-1-1-0-1-0 1450, where the first 7 num- 
bers characterize the metronome bobs position relative to 
the photo-cell (whether the gate is open or closed) and 
the ninth number is the time, where one time unit cor- 
responds to 64 microseconds. Since we are interested in 
the dynamics of this system from the viewpoint of syn- 
chronization, for numerical evaluation we will compute 
the classical order parameter, r, of the Kuramoto model 
i: 

r exp(i0) = ^ X! Gxp(i6lj). (1) 

j 

Here is an average phase of the whole ensemble, 6j is 
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FIG. 3: (a) The experimental setup, with the 
metronomes placed on the platform and the wiring that 
carries information on the position of the metronomes ' 
phases, (b) One of the used light-gates (Kinghright 
KTIR 0611 S), composed of an infrared LED and a 
photo-transistor. 



the phase of the j-th metronome, A'' is the number of 
metronomes, and i is the imaginary unit. 

From the recorded data we can get directly only those 
moments in which the metronomes' bobs pass through 
the light-gates, so some additional steps arc needed in or- 
der to get the 6j phases of all metronomes and to compute 
for given time moments the Kuramoto order-parameter. 
In order to achieve this, first we excluded from the data 
those time-moments when the metronome's bob passed 
the second time through the light-gate in one period, and 
after that we remained with the pass-times correspond- 
ing to a given directional motion only. With this " cleaned 
data" we calculated the period of each cycle and interpo- 
lated this time-interval for the 9j phases (between and 
27r, corresponding to the state of a Kuramoto rotator), 
assuming a uniform angular- velocity. In such manner in 
each time-moment, the Oj phase of each metronome (con- 
sidered here as a rotator) can be uniquely determined and 
the Kuramoto order parameter ([T|) was computed. 



uj{BPM) 


160 


168 


176 


184 


192 


208 


(7 {BPM ■ IQ-'') 


8.4 


7.9 


7.8 


9.8 


8.5 


8.7 



TABLE I: Standard deviation of the used seven 
metronomes for different frequencies 



Before starting the experiments we monitored sepa- 
rately each metronome, and recorded their exact fre- 
quency for all the standardly marked rhythms. In such 
manner we have selected those 7 metronomes that had 
their standard frequencies relatively close to each other, 
and precisely measured these values. From these values 
the standard deviation of the metronomes' natural fre- 
quencies could be determined (Table |l| . 

III. EXPERIMENTAL RESULTS 

As already described in the introductory section, the 
metronomes oscillate with different natural frequencies, 
depending on the position of the adjustable weight on the 
metronomes' rod. For our experiments we have used the 
standard frequencies marked on the metronome. These 
frequencies are given in BPM units. 

Before discussing in detail the experimental results, 
we have to emphasize, that independently of the cho- 
sen initial condition, only in-phase synchronization of the 
metronomes was observed. Reasons for this will be given 
in a separate section (Section |Vl|. 

In the very first experiments we were interested 
in studying how the chosen frequency influences the 
detected synchronization level. We fixed all the 
metronomes' frequencies on the same marked value and 
placed them symmetrically on the perimeter of the ro- 
tating platform as indicated in Figure [3^. Naturally, this 
does not mean that in reality their frequencies were ex- 
actly the same, since no two real macroscopic physical 
systems can be exactly identical. We initialized the sys- 
tem by starting the metronomes randomly, and let the 
system composed by the metronomes and platform freely 
evolve. 

For each considered frequency value we made 10 mea- 
surements, collecting data for 10 minutes. The dynamics 
of the computed Kuramoto order parameter averaged on 
the 10 independent realizations are presented on Figure 

Hi. 

The results suggest that by increasing the metronomes' 
natural frequencies, the degree of the obtained synchro- 
nization level increases. The standard deviations of the 
natural frequencies of the independent oscillators is indi- 
cated in Table |T1 

Since there is no clear trend in this data as a function 
of cj, the obtained result suggests that the observed effect 
is not due to a decreasing trend of the metronomes' stan- 
dard deviation. We have also found that for the standard 
metronome frequencies below 160 BPM the system did 
not synchronize. It is interesting to note however, that if 



4 




° co= 


:160BPM 


° co= 


:168BPM 


» co= 


:176BPM 


" co= 


:184BPM 


+ co= 


:192BPM 


X co= 


:208BPM 




0.3, 
b) 







40 150, 160 170 

.1 , , t(s), ^ 

m 



"TOD" 

t(s) 



300 



FIG. 4: Dynamics of the order parameter for different 
natural frequencies and metronome numbers, (a) 
Results for seven metronomes, different curves 
corresponding to different frequencies as indicated in the 
legends, (b ) Results are for the same natural frequency 

(^192 BPM) and different metronome number as 
indicated in the legend. On both graphs the results are 
averaged on 10 independent measurements. 
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3 


4 


5 


6 


7 


(7 [BPM ■ IQ-'^) 


5.1 


8.1 


7.5 


7.1 


6.7 


8.5 



TABLE II: Standard deviation of the metronomes 
natural frequencies for different metronome number on 
the rotating platform. 



get for the detected synchronization level a clear trend: 
increasing the number of metronomes will result in a de- 
crease of the synchronization level. 



IV. THEORETICAL MODEL 

Inspired by the model described in (T^], it is possible 
to consider a simple mechanical model for the system in- 
vestigated here. The model is composed of a rotating 
platform and physical pendulums attached to its perime- 
ter, as it is sketched on Figure [51 




one inspects only visually or auditory the system, he/she 
would observe no synchronization for frequencies already 
below 184 BPM. This means that we are not suited to 
detect partial synchronization with an order parameter 
below r = 0.75. 

In a second experiment we were investigating the in- 
fluence of the metronome number on the synchroniza- 
tion level. In order to study this, we have fixed the 
metronomes on the same frequency (192 BPM) and re- 
peated the previous experiment with increasing number 
of metronomes placed on the rotating platform. 

Again, for obtaining accurate results we performed 10 
measurements for each configuration, and averaged the 
observed order parameter. The averaged results are pre- 
sented on figure [3)3. 

Although the standard deviation of the metronomes' 
natural frequencies (Table [TT| does not present a clear 
trend as a function of the metronome number, iV, we 



FIG. 5: Schematic view and notations for the considered 
mechanical model. The white dot denotes the center of 
mass of the physical pendulums and the gray dot is the 
suspension axis. 
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The Lagrange function of such a system writes as 
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\d_ 

di' 



2 ^ 

(jii cos I — mighjil — cos Oj 



+ 



(2) 



i=l 



The first term is the kinetic energy of the platform, the 
second is the kinetic energy due to the rotation of the 
pendulum around its center of mass, the third one is 
the kinetic energy of the pendulum's center of mass, and 
the last term is the potential energy of the pendulum. 
In the Lagrangian we have used the following notations: 
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the index i denotes the pendulums, J is the moment of 
inertia of the platform with the metronomes on it - taken 
relative to the vertical rotation axes, is the angular 
displacement of the platform, rrii is the total mass of the 



pendulum {im « W^^ + neglecting the mass of the 
rod), hi is the distance between the center of mass and 
the suspension point of the pendulum, xi is the horizontal 
displacement of the center of mass of the pendulums due 
to the rotation of the platform, 6'j is the displacement of 
the i-th pendulum's center of mass, in radians, Ji is the 
moment of inertia of the pendulum relative to its center 
of mass and Wi is the angular velocity of the rotation of 
the pendulum relative to its center of mass. It is easy to 
realize that Xi = R4> and uji = 9i. Assuming now that the 
mass of all the weights suspended on the metronomes bob 
are the same {w[^'^ = wi, Wj*'' = W2, and consequently 
rrii = m), and disregarding the niighi constant terms, 
one obtains: 



L' = 



J + NmR^ 



iRcf) hi 



E 



cos 9i ■ 6i 



mg hi cos 6i (3) 



The Euler-Lagrange equations of motion yields: 



(J + NmR^)ij) + mR^h^ [0, 



■ sin di 



(4) 



\mhl + Ji\Bi + mR(j)hi cos 9i + mghi sin 9i = 0. (5) 

The above equations of motion are for a Hamiltonian 
system where there is no damping (no friction) and no 
driving (excitation). Friction and excitation from the 
metronomes' driving mechanism has to be taken into ac- 
count with some extra terms. The system of equations 
of motion writes as 

(J + NmR^)ij) + mRY] h, [Oi cos 9^ - 9^ sin 9,] + 



+C00 + ^M, =0 



(6) 



also that the excitation is given in the good direction 
(direction of the motion). It is easy to see that the total 
momentum transferred, Wltransi to the metronomes in a 
half period (T/2) is always M: 



t+T/2 



MS{9i)d9, = M 




This driving will be implemented in the numerical solu- 
tion as 



if 9,{t -dt) <0 and 9i{t) > 
if 9,{t-dt) > and 9i{t) < 
in any other case 



where dt is the time-step in the numerical integration of 
the equations of motion. It is immediate to realize that 
this driving leads to the same total momentum transfer 
M, as the one defined by equation 

The coupled system of equations (j6l7|) can be written 
in a form more suitable for numerical integration: 



iRJ2i f^i^i sin^i ^ - J2i 



A + B-C 



D 



mR(f>hi cos 9i — mghi sin 9i — cg9i 
mhf + Ji 



(9) 



(10) 



where 



A 2 n sin 9i cos 6*, 

A = m gR ^ ■ 



B = mRce 



mh'f + Ji 
X - hi9i cos9i 



C = mR 



^ mhf + Ji ' 
y~v hiMi COS 9i 



m 



hf + jr 



D = {J + NmR^ 



hf cos^ 9i 
mhf + Jj 



[mhf + J,]9, 



mR(j)hi cos 9i 
+mghi sin 9i 



(7) 



where and eg are coefficients characterizing the friction 
in the rotation of the platform and pendulums, and 
are some instantaneous excitation terms defined as 



MSi9,)9,, 



(8) 



where i5 denotes the Dirac function and M is a fixed 
parameter characterizing the driving mechanism of the 
metronomes. The choice of the form for M, in equa- 
tion ([5]) means that excitations are given only when the 
metronome's bob passes the 9 — Q position. The term 9 
is needed in order to ensure a constant momentum input, 
independently of the metronomes amplitude. It insures 



Now taking into account that the metronomes' bobs have 
the form sketched in Figure [2]3 and the Li distances are 
fixed and assumed to be identical for all the metronomes, 
the hi and J,; terms of the physical pendulums in our 
model will be calculated as: 



1 



J, = wiLl+W2{LfY. 



Wi + W2 



(11) 

(12) 



V. REALISTIC MODEL PARAMETERS 

The parameters were chosen in agreement with the 
characteristic values of our experimental device: wi = 



0.025 kg, W2 = 0.0069 kg, Li = 0.0358 m, La £ 
[0.019, 0.049]m depending on the chosen natural fre- 
quency, R = 0.27 m and J G [0.0729,0.25515]% 
depending on the number of metronomes placed on the 
platform. 

The damping and excitation coefficients were esti- 
mated as follows. For the estimation oi cg, a single 
metronome on a rigid support was considered. Switching 
off the excitation mechanism, a quasi-harmonic damped 
oscillation of the metronome took place. The exponen- 
tial decay of its amplitude as a function of time uniquely 
defines the damping coefficient, hence a simple fit of the 
amplitude as a function of time allowed the determina- 
tion of eg. Switching now the exciting mechanism on, 
leads the metronome to a steady-state oscillation regime 
with constant amplitude. Since cg has already been mea- 
sured, this amplitude value is defined by the excitation 
coefficient M. Solving equations © and (I7|) for a single 
metronome and tuning it until the same steady-state am- 
plitude is obtained as in experiments makes possible the 
estimation of M. Now that both cg and M are known, the 
following scenario is considered: all the metronomes are 
placed on the platform and synchronization is reached. 
Then the platform has a constant-amplitude oscillatory 
motion. In order to determine c^, its value is tuned while 
solving the equations ©-(I?]) until the same amplitude 
of the disc's oscillations is obtained as in experiments. 
This way all the parameters in the model can be re- 
lated to the experimental quantities. Using the above 
defined method, we estimated the following parameter 
values: cg = 0.00005 kg m^/s, = 0.00001 kg m'^/s 
and M = 0.000095 Nm/s. 
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FIG. 6: Dynamics of two metronomes as a function of 
time, and the quantities used for defining the z order 
parameter 



appropriate for such purpose, since it cannot distinguish 
between the totally unsynchronized and the anti-phase 
synchronized dynamics. In order to introduce a proper 
order parameter, let us consider the dynamics of two 
metronomes as a function of time by plotting 0i.2(t) (Fig- 
urc[6]). Let us denote by if"" and t^"'^ the time-moments 
where metronome 1 has local minimum and maximum 
6*1 values, respectively. We denote by T™"^ the time- 
moments where metronome 2 has local maximum 82 val- 
ues. With these notations, we define the following two 
time-like quantities, that are characterizing the average 
time- interval of the maximum position of the 02{t) rel- 
ative to the maximum and minimum positions of 9i{t), 
respectively: 



VI. IN-PHASE SYNCHRONIZATION VERSUS 
ANTI-PHASE SYNCHRONIZATION OF TWO 
METRONOMES 

As it is described in the introductory section, in the 
case of two coupled oscillators many previous works have 
reported a stable anti-synchronized state [l^, [ill ■ Due 
to the fact that independently of the starting condition 
no such stable phase was observed in our experiments, we 
feel that investigating this issue is important. Starting 
from our theoretical model described in section IIVI we 
will show that the in-phase synchronization is favored, 
whenever there are equilibrated damping and driving 
forces acting on the metronomes. 

Let us investigate first the case without any damp- 
ing and with no driving forces. The equations of mo- 
tion for such a system are given by equations ((4]) and 
([5]). Considering the case of two metronomes (A^ = 2) 
with small-angle deviations {0^2^ « 7r/2), we inves- 
tigate the synchronization properties of such a system. 
The synchronization level will be studied here by an ap- 
propriately chosen order parameter for two metronomes, 
z, that indicates whether we have in-phase or anti-phase 
synchronization. The Kuramoto order parameter is not 



ti = (min{|tr 
t2 = {mm{\t 



rpmax 
j 



min rpmax 

~ j 



(13) 
(14) 



In the equations from above the averages are consid- 
ered for all j maximum positions of 62 (t) , and the " min" 
notation refers to the minimal value of the quantity in 
the brackets. Now, the z order parameter is defined as: 



t2-h 

t2+tl 



(15) 



It is easy to realize that z is bounded between —1 and 
1 . For a totally in-phase synchronized dynamics we have 
ti = 0, leading to 2; = 1. For a totally anti-phase syn- 
chronized dynamics t2 = 0, and we get z ~ —\. Negative 
z values suggest a dynamics where the anti-phase syn- 
chronized states are dominant, positive z values suggest 
a dynamics with more pronounced in-phase synchronized 
states. 

The z order parameter was estimated numerically for 
different initial conditions. A velocity Verlet-typc algo- 
rithm was used, and simulations were performed up to a 
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t = 10000 s time interval, with a dt = O.l s time-step. 
(We remark, that this time-step is sufficiently small, since 
for dt ~ 0.001 s we obtained very similar results.) 

Initially the deviation angle of the first metronome was 
chosen as 9i{0) = 9 max =0.1 rad and (0) was chosen 
in the interval [—0.1,0.1] rad, leading to various initial 
phase- differences between them. The obtained z values 
as □ 




- J=0.01 
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FIG. 7; The z order parameter as function of the initial 
phase ^2(0) of metronome 2. (a) Simulation results 
without damping (cg = 0, c^j, ~ 0) and no driving 
(M — Q). (h) Simulation results with the realistic 
damping and driving parameters given in Section [2 

The obtained results suggest that for the friction-free 
case with no driving (Figure [7^), depending on how the 
dynamics is initialized, both the synchronized and anti- 
synchronized states are observable. If the anti-phase ini- 
tial condition is closer, the dynamics will be mostly anti- 
phase synchronized. Consequently, if the initial condi- 
tion is closer to the in-phase synchronization, the dy- 
namics will be mostly in-phase synchronized. Increasing 
the inertial momenta of the platform will lead to a better 
synchronization level. In case both damping and driving 
is present and the parameters are the ones experimen- 
tally measured in the system (Figure [7|d), the anti-phase 
synchronization is suppressed and only in-phase synchro- 
nization appears. Anti-phase synchronization is unlikely 



to be observed, this will appear only in the case when 
the two metronomes are started exactly in anti-phases 
(02(0) ~ -0i(O)). 

In view of these results, one can understand why in 
our experiments only the in-phase synchronized dynam- 
ics was observed. 



VII. SIMULATION RESULTS FOR SEVERAL 
METRONOMES 

The equations of motion (|9ll0p were numerically in- 
tegrated using a velocity Verlet-type algorithm as inte- 
gration method. A time-step oi dt — Q.\ s was chosen. 
As a test of the model, first we tried to reproduce the 
experimental results. Our first simulation is intended to 
reproduce the experimental results presented in Figure |4l 
Seven metronomes with the same natural frequencies as 
the experimentally measured ones were considered, and 
the time-evolution of the Kuramoto order parameter was 
computed. Results for different frequency values are pre- 
sented in Figure [5^1. 

Except for the first maximum, the trends are in good 
agreement with the experimental results from Figure 2] 
We believe that this maximum appears due to the differ- 
ences in the initialization of the dynamics in the experi- 
ments and simulations. 

As a second simulation experiment, we have studied 
the time-evolution of the order parameter for different 
pendulum number, setting the same 192 BPM natural 
frequency as in the experiments. Again, similarly to the 
experiments, we averaged the results for 10 independent 
simulations. The obtained trend is sketched on Figure 

The trend of the simulation results is in agreement 
with the experimental ones: increasing the metronomes' 
number results in a decrease in the observed synchroniza- 
tion level. In simulations however, this decrease is not as 
evident as in the experiments. We believe that the dif- 
ference is due to the simplified manner how the slight 
differences between the metronomes are handled in the 
simulations. We recall that in the model, the differences 
between metronomes is included only via the Lj*^ terms. 
In reality the other parameters of the metronomes can 
also differ. 

These simulation results suggest that we have a rea- 
sonable model for understanding the synchronization in 
such a system. One can investigate thus by simulations 
several interesting cases, that are not feasible in the ex- 
periments. Many interesting questions can be formulated 
in such manner. In the following we will focus however 
only on clarifying the problems that we have experimen- 
tally investigated, namely the influence of the number of 
oscillators and the chosen natural frequency on the ob- 
served synchronization level. Computer simulations will 
allow us to consider a high number of oscillators and will 
also allow a continuous variation of the metronomes' nat- 
ural frequencies. Particularly, we are interested in clari- 
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FIG. 8: Simulation results for the time-evolution of the 
Kuramoto order-parameter, (a) The curves represent 
results for the same natural frequencies, as the ones 
used in experiments, (b ) Results for the same pendulum 
number as the ones used in experiments. For both 
graphs the presented results are an average over 10 
independent simulations. 



fying how the stable synchronization level is varying as a 
function of the metronomes' number, and how the char- 
acteristic time needed to achieve synchronization varies 
as a function of the metronomes' number. 

First we have reconsidered the results presented in Fig- 
ure [HId and investigated the system for large number of 
metronomes. Also, for the sake of a better accuracy we 
considered an average over 100 simulations for each fixed 
values. The same 192 BPM natural frequency was 
fixed for all oscillators. Simulations with up to = 256 
metronomes were considered and the obtained results are 
presented on Figure [Hh. 

The results confirm our previous conclusion: increas- 
ing the metronomes number will result in a decrease in 
the observed synchronization level. However, a quantita- 
tive analysis of the results suggests that in the limit of 
N ^ oo (we considered up to 1024 metronomes) the Ku- 
ramoto order parameter converges to a non-zero value. 
Comparing the results from Figure IH^ with the ones pre- 
sented in Figure [B)d confirms also what one would natu- 
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FIG. 9: Simulation results for the equilibrium 
synchronization level, (a) Time- evolution of the r order 
parameter for different pendulum number, as indicated 
in the legend. One graph is obtained after an average on 
100 independent simulations, (b) The r order parameter 

as a function of the pendulum number. The circles 
indicate the simulation results and the continuous curve 
shows the fit 116\) with parameters: a — 0.041, b — 0.049 
and c ^ 0.971. 



rally expect based on the Kuramoto model: considering 
totally identical metronomes with no standard deviation 
of their natural frequencies, leads to an increased Ku- 
ramoto order parameter. Simulations suggest that the 
values of the stabilized order parameter converge expo- 
nentially to a stable c = 0.971 value 

r(Ar) = aexp(-feA') -t-c, (16) 

with the best fit parameters a = 0.041, b — 0.049. On 
Figure [9)d we plot this A^ dependence, and indicate also 
the fit given with the above mentioned parameters. 

The results from above prove that by increasing the 
number of metronomes placed on the rotating platform 
the partial synchronization survives. The system con- 
verges to a dynamic state with a non-zero synchroniza- 
tion order parameter. We have also learned from Fig- 



ure [9] that the equihbrium value of the order parame- 
ter does not vary too much by increasing the number of 
metronomes. In contrast with this small variation, the 
time needed to achieve synchronization increases more 
evidently. In order to define a characteristic time needed 
to achieve synchronization, we define an averaged order 
parameter, q, that is suitable to capture the initial part 
of the time-evolution: 

9=/*"°'^^*, (17) 

Jdt ' 




Here dt is the integration time-step for the equations 
of motion and tmax is roughly the time-length of the con- 
sidered simulation. 

This order parameter is a weighted average on the 
[dt,tmax] interval, where the weight factors are favor- 
ing the initial time-moments. As a result of this, this 
order parameter captures the initial transient part of 
the dynamics. When the value of the Kuramoto order 
parameter does not change significantly (the case when 
the number of metronomes, N, were varied) the q order 
parameter is appropriate for revealing the dynamics in 
which synchronization is achieved quicker. Consequently, 
a characteristic time, tc, to achieve synchronization can 
be defined as: = l/q. Naturally, the value of q and tc 
depend on both the dt and tmax values (the limits in the 
integral pT)) ). The good news however is, that although 
the values are dependent on the chosen dt and tmax val- 
ues, seemingly there is a scaling of g as a function of 
N: 



C ■ N~^'^ + K{tmax) 



(18) 



Results confirming this scaling are presented in Figure 
[TOl Simulations yield the C = 5.55 constant. Experimen- 
tal results plotted on Figure ITOb indicate the same scal- 
ing. Moreover, the slope of the linear fit gives C = 5.38 
in excellent agreement with the simulation results. 

Computer simulations allow us also to investigate in 
more detail the dependence of the r and q order pa- 
rameter as a function of the natural frequencies of the 
pendulums. Results for different pendulum number are 
given in Figure [TT] These results confirm again, that by 
decreasing the metronomes' natural frequency the syn- 
chronization level decreases. Increasing the metronome 
number and/or by decreasing their natural frequencies 
makes harder and harder to achieve a synchronized state. 



VIII. CONCLUSIONS 

The dynamics of a system composed of coupled 
metronomes was investigated both by simple experiments 
and computer simulations. We were interested in finding 
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FIG. 10: Scaling of the q order-parameter, (a) Symbols 
are simulation results for various t^ax values, and the 
curves are linear fits. One can observe that the slopes 
are roughly the same, and the linear fit is a good 

approximation, (b ) Scaling for the experimental results. 



the conditions for an emerging synchronization. Contrar- 
ily with many previous studies, here the problem was an- 
alyzed not from the viewpoint of dynamical systems, but 
from the viewpoint of collective behavior and emerging 
synchronization. 

The experiments suggested that there is a limiting nat- 
ural frequency of the metronomes below which sponta- 
neous synchronization is not possible. By increasing the 
frequency above this limit partial synchronization will 
emerge. The obtained synchronization level increases 
monotonically with increasing the natural frequency of 
the oscillators. The experiments proved that by increas- 
ing the number of metronomes in the system the final 
stable synchronization level drops and the time needed 
to achieve the final dynamical equilibrium state increases. 

In order to better understand the dynamics of the sys- 
tem a realistic model was build. The parameters of the 
model were fixed in agreement with the experimental 
conditions and the equations of motion were numerically 
integrated. The model proved to describe well the ex- 
periments, and reproduced the experimentally observed 
trends. The model allowed a fine verification of our find- 
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FIG. 11: Simulation results for the influence of the 
metronomes natural frequencies lo on the obtained 
synchronization level, measured by the r (Figure a) and 
q (Figure b) order parameters. Different curves are for 
different metronome number, as indicated in the legends. 



ings regarding the conditions under which spontaneous 
synchronization emerges. It was found that by increas- 
ing the number of metronomes in the system the synchro- 
nization level converges exponentially to a stable limit. 
Also, it was found that the time needed to achieve syn- 
chronization scales in a universal manner with the num- 
ber of metronomes placed on the rotating platform. 

The successes of the elaborated model opens an easy 
way for many further studies regarding the dynamics 
of this system. Indeed, many other interesting ques- 
tions can be formulated regarding the influence of the 
metronomes and rotating platform parameters on the 
obtained synchronization level and the observed trends. 
Also, one can study systems where the metronomes or 
groups of metronomes are fixed to different natural fre- 
quencies, or where there is an external driving on the sys- 
tem. The elaborated model has the advantage that the 
equations of motion are easily integrable and the model 
parameters are realistic having a direct connection to the 
parameters of an experimentally realizable system. 

Finally, we hope that the novel experimental setup and 
the results presented here will help in clarifying some as- 
pects for one of the oldest problem in physics, namely 
the spontaneous synchronization of coupled pendulum 
clocks. Although several similar problems have been con- 
sidered in previous studies, we have shown that there are 
still many fascinating aspects that one can investigate in 
such simple mechanical systems. 
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